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Abstract — Mutually unbiased bases (MUBs) are a primitive 
used in quantum information processing to capture the principle 
of complementarity. While constructions of maximal sets of d + 1 
such bases are known for system of prime power dimension d, it is 
unknown whether this bound can be achieved for any non-prime 
power dimension. In this paper we demonstrate that maximal sets 
of MUBs come with a rich combinatorial structure by showing 
that they actually are the same objects as the complex projective 
2-designs with angle set {0, 1/d}. We also give a new and 
simple proof that symmetric informationally complete POVMs 
are complex projective 2-designs with angle set {l/(d+l)}. 

I. Introduction 

Two quantum mechanical observables are called comple- 
mentary if and only if precise knowledge of one of them 
implies that all possible outcomes are equally probable when 
measuring the other, see for example [19, p. 561]. The 
principle of complementarity was introduced by Bohr [6] in 
1928, and it had a profound impact on the further development 
of quantum mechanics. A recent application is the quantum 
key exchange protocol by Bennett and Brassard [3] that 
exploits complementarity to secure the key exchange against 
eavesdropping. 

We mention a simple mathematical consequence of this 
complementarity principle, which motivates some key notion. 
Suppose that O and O' are two hermitian d x d matrices 
representing a pair of complementary observables. We assume 
that the eigenvalues of both matrices are multiplicity free. 
It follows that the observables O and O' respectively have 
orthonormal eigenbases B and B' with basis vectors uniquely 
determined up to a scalar factor. 

The complementarity of O and O 1 implies that if a quantum 
system is prepared in an eigenstate b' of the observable O', 
and O is subsequently measured, then the probability to find 
the system after the measurement in the state b € B is given 
by |(fe|6')| 2 = 1/d. Recall that two orthonormal bases B and 
B' of C d are said to be mutually unbiased precisely when 
|(6|6')| 2 = l/d holds for all b e B and b' e B' . Thus 
the eigenbases of non-degenerate complementary observables 
are mutually unbiased. Conversely, we can associate to a 
pair of mutually unbiased bases a pair of non-degenerate 
complementary observables. 

There is a fundamental property of mutually unbiased bases 
that is invaluable in quantum information processing. Suppose 
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that we want to determine the density matrix p of an ensemble 
of quantum systems using as few non-degenerate observables 
as possible. We assume that it is possible to make a com- 
plete measurement of each observable O = J2beB Xb l> 
meaning that the statistics tr(p \b) (b\) = (b\ p \b) is known for 
each eigenvalue Xb in the spectral decomposition. Ivanovic 
showed in [12] that complete measurements of at least d + 1 
observables are needed to reconstruct the density matrix. He 
also showed that this lower bound is attained when d + 1 non- 
degenerate pairwise complementary observables are used. 

A simple example is provided by the Pauli spin matrices a x , 
a y , a z . A complete measurement of these three observables 
allows to reconstruct a 2 x 2 density matrix, a fact appar- 
ently known to Schwinger [18]. Nowadays, we know how to 
do this state tomography process — at least in principle — in 
dimensions d = 3,4, and 5. It is an open problem whether 
it is possible to perform this kind of state tomography in 
dimension 6, because the construction of a set of 7 mutually 
unbiased bases in dimension d = 6 is elusive. 

II. Mutually Unbiased Bases 

Definition 1: Two orthonormal bases B and C of C d are 
called mutually unbiased iff |(&|c)| 2 = 1/d holds for all b € B 
and ceC. 

The goal is to construct d+1 mutually unbiased bases 
(MUBs) in any dimension d > 2. There are several con- 
structions known to obtain MUBs. At least for prime power 
dimension the problem is completely solved. This follows 
from Constructions I-III below. However, in any dimension 
other than a prime power it is unknown if a maximal set 
of d + 1 MUBs can be found. The best known result is 
Construction IV below which only works in dimensions d 
which are squares and never gives a maximal set of MUBs. 

Construction I (Wootters and Fields [24]) Let q be an odd 
prime power. Define 

|„, \ „— 1/2 /, ±r(ax 2 +bx)\ ~- /~iq 

\Va,b) — q [u p K ')ieF, e L', 

with uj p = exp(27ri/p). Then the standard basis together with 
the bases B a = {|v a ,h) \b G F f/ }, a € F f/ , form a set of q + 1 
mutually unbiased bases of C q . 



Construction II (Galois Rings [13]) Let GR(4, n) be a finite 
Galois ring with Teichmiiller set T n . Define 

K, b ) = 2-"/ 2 (cxp ^tr(a + 2b)^j ) ^ . 

Then the standard basis together with the bases M a = 
{\ v a,b) \b E T n }, a E T n , form a set of 2" + 1 mutually 
unbiased bases of C 2 . 

Construction III (Bandyopadhyay et al. [1]) Suppose there 
exist subsets C\ , . . . , C m of a unitary error basis B such that 
\Ci\ = d, d n Cj = {Id} for i ^ j, and the elements of d 
pairwise commute. Let Ms be a matrix which diagonalizes Ci. 
Then Mi,, ... , M m are MUBs. 

Construction IV (Wocjan and Beth [23]) Suppose there are 
w mutually orthogonal Latin squares [4], each of size d x d 
over the symbol set S = {1, ...,d}. Then w + 2 MUBs 
in dimension d 2 can be constructed as follows. With each 
Latin square L (and additionally the square (1, 2, . . . , n)* <8> 
(1,...,1)) we can associate vectors of length d over the 
alphabet {1, . . . , g? 2 }: for each symbol a E S define a vector 
sl,o € C d as follows: start with the empty list SL, a = 0- 
Then traverse the elements of L column-wise starting at the 
upper left corner. Whenever a occurs in position in L, 
then append the number i + jd to the list sl.o.- The other 
ingredient to construct these MUBs is an arbitrary complex 
Hadamard matrix H = (hij) of size d x d. For each Latin 
square L and each a,j E {1, ...,d} define a normalized 
vector \vL,a,j) ■= e s L a [i]hi,j> where e, are the 

elementary basis vectors in C d . Then the bases given by 
Bl '■= {\vL.aj) ■ ct,j = 1, . . . , d} together with the identity 
matrix 1^2 form a set of w + 2 MUBs. 

Example 1: In dimension d = 3 Construction I yields the 
bases 

3-V2{ (1,1,1), (1,0,3, U,f), (1,W§,«3)}, 

3- 1 /a{(i jW3)W3 ) ) (1,^ 2 ,1), (1,1, c^)}, 

3- 1 /2{(l, w 2, w 2), (1,0,3,1), (1,1, W 3 )}, 

which together with the standard basis I3 form a maximal 
system of four MUBs in C 3 . 

Example 2: In dimension d = 4 Construction II yields the 
bases (where we have abbreviated " + " for 1 and " — " for 
— 1 and i = 

!{(+,+,+,+),(+,+,--),(+,--+),(+,-+,-)}, 

|{(+, -, ~i, -i), (+, -, i, i), (+, +, -«)> (+) +> ~h 

-*i -i, -)i (+, -ij *i +), (+, ij *i -), (+, i, ~i, +)}, 

|{(+, -ii -, -«)) (+> +> *)) (+> +> _i )> (+> -> *)}• 
These four bases and the standard basis I4 form an extremal 
set of five MUBs in C 4 . 

A basic question is how many bases can be achieved in 
general dimension. To this end, we define the function M : 
N — * N as follows: 

M(n) := max{|S| : B is a set of MUBs in C"} 



Then we have that: 

• M(p r ) — p r + 1 for p prime, r E N, 

• M (ri) < n + 1 for all n e N, 

« M(mn) > min{M(m), M(n)} for all m, n E N. 

• M(d 2 ) > N(d), where AT(d) is the number of mutually 
orthogonal Latin squares of size d x d. 

An open problem is to show that liminfn^oo M(n) = 00. 

III. Welch's Lower Bounds 

Suppose that X is a finite nonempty set of vectors of unit 
norm in the complex vector space C d . The vectors in X satisfy 
the inequalities 

iij2 E Kz|y>r>73+iUy (D 

1 1 1,1/ex v k 1 

for all integers k > 0. Welch derived these bounds in [22] 
to obtain a lower bound on the maximal cross-correlation of 
spreading sequences of synchronous code-division multiple- 
access systems. Blichfeld [5] and Sidelnikov [20] derived 
similar bounds for real vectors of unit norm. 

A set X attaining the Welch bound Q for k = 1 is 
called a WBE-sequence set, a notion popularized by Massey 
and Mittelholzer [15] and others. Using equation Q, it is 
straightforward to check that the union of d + 1 mutually 
unbiased bases of C d form a WBE-sequence set. These 
extremal sets of mutually unbiased bases are even better, since 
they also attain the Welch bound for k = 2. In fact, we 
show that a sequence set attains the Welch bounds Q for all 
k < t if and only if it is a i-design in the complex projective 
space CP d -\ 

Let us introduce some notation. Let S d ~ 1 denote the sphere 
of unit vectors in the complex vector space C d . We say that 
two vectors u and v of S^ 1 are equivalent, in signs u = 
v, if and only if u = e l6 v for some € R. It is easy to 
see that = is an equivalence relation. We denote the quotient 
manifold S"^ 1 /^ by CS' 1 - 1 . Notice that the manifold CS^ 1 
is isomorphic to the complex projective space CP d_1 , but we 
prefer the former notation because normalizing vectors to unit 
length is common practice in quantum computing. 

Lemma 1: Let p be the unique normalized J7(d)-invariant 
Haar measure on the complex sphere CS''' -1 . For any x E 
we have 

/ My)\ 2k dp(y) = - [ ^ r 

Jcs*-i ( k ) 

Proof: The unitary group U[d) acts transitively on the 
manifold C5 d_1 . This means that for any y E CS d ^ 1 there 
exists a unitary matrix U mapping y to the first basis vector, 
Uy = e\. Therefore, 

f \(x\y)\* k d»(x) = [ \(Ux\ ei )\ 2k dn(x) 
Jcs d - 1 JCS 11 - 1 

= / Kx|ei)| 2 *dA*(x), 
Jcs d - 1 

where the last equality holds because of the [/(d)-invariance 
of the measure p. Using Proposition 1.4.9 from Rudin [17], 



we obtain 

which proves the claim. 



cs^- 1 



x k \ 2 dfi(x) 



IV. Complex Projective i-DESiGNS 

We now present some background material on complex 
projective designs. We will relate those later on to the systems 
of vectors formed by a maximal set of MUBs. 

Let us first introduce some notation. We denote 
by Hom(fc,£) the subset of the polynomial ring 
C[xi, ... ,x<i,yi, ■■■ ,Vd] that consists of all polynomials that 
are homogeneous of degree k in the variables xi,...,Xd 
and homogeneous of degree £ in the variables yi, ...,j/<j. 
We associate to each polynomial p in Hom(fc,£) a function 
p on the sphere 5 d_1 by defining p (Q = p(£, £) for 
£ G S d ~ 1 . The function p Q is called the "restriction" of p 
onto the complex sphere. It follows from the homogeneity 
conditions of the polynomial p that p (e vi> £ i ) = e ( k ~ ( >p (£) 
holds for all d G R. Therefore, we obtain a well-defined 
polynomial function on CS' d_1 only if k = £. We define 
Hom(fc, k) Q = {p : p G Hom(/c, &)}. 

Definition 2: A finite nonempty subset X of C5 d_1 is a 
i-design in CS d ~ 1 iff the cubature formula 



n^E^ 
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M(CS* d - 



cs d - 1 



f(x)dn(x) 



xex 

holds for all / in Hom(i,i) . 

We now show a characterization of i-designs in terms of 
the inequalities by Welch given in equation Q. 

Theorem 1: Suppose that X is a finite nonempty subset of 
CS d_1 . Then the following statements are equivalent: 

1) The set X is a £-design in OS"* -1 ; 

2) for all x G C d and all k in the range < k < t, we have 
the equality 

(x\x) k 



rd+k-l\ \ X \ ^ 



\{x\y)\ 



2k. 



(2) 



3) the set X satisfies the Welch bounds Q with equality for 
all k in the range < k < t, that is 

1 ^ ou 1 



\x\ 



E \^y)\ 2k 



0<k<t. (3) 



i+k-l\ ' 
x,y£X V k I 

Proof: We show that 1) implies 2). Fix a vector x € C d . 
Note that p(y) = |(a;|y)| 2fe = (y\x) k (x\y) k is a polynomial 
function in Hom(fc,fc) . Since X is a t-design, the exact 
cubature formula 



m£'<*>' 2 *=L 



\(x\y}\ 2k d^y) 



holds for all k in the range < k < t. By Lemma^ the latter 
integral evaluates to ( c 
(|3 holds for all k <t. 



integral evaluates to ( d+k . X ) , which proves that equation 



We show next that 2) implies 3). We observe that Q holds 
for all k < t, hence summing over x G X yields (|3j- 

Finally, we show that 3) implies 1). Suppose that equation 
(0 holds. For a vector x G C d , we denote by x® k the k- 



fold tensor product x 



x (8) 



X G C Q 



Note that 



Consider the d -dimensional vector 



\x I E 



x® k ®x® k dn(x). 



CS d 



Evaluating the inner product of £ with itself yields 



LW E K a 



2 k 



x,y£X 



Kxly)!^^^)^^), (4) 



which is equal to > 0. The inner integral evaluates to 

( d+ t _1 ) Lemma[2 and the double integral has the same 
value, because the measure fi is normalized. It follows from 
our assumption that the right hand side vanishes. By con- 
struction of £, we can conclude that averaging over X yields 
an exact cubature formula for all monomials in Hom(fc, fc) , 
hence, by linearity, for all polynomials in Hom(fc, k) - This 
means that X is a i-design. ■ 

Remark 1: Equation provides a short proof of the Welch 
inequalities Q. The analogue for real spherical i-designs 
of the above result is sketched in [7]. A connection to the 
existence of certain isometric Banach space embeddings is 
given in [14]. 

V. Uniform Tight Frames 

A finite subset F of nonzero vectors of C d is called a frame 
if there exist nonzero real constants A and B such that 

A\\ir<J2Kf\v)\ 2 <B\\v\\ 2 
feF 

holds for all v G C d . The notion of a frame generalizes the 
concept of an orthonormal basis. The linear span of the vectors 
in F the space C d , but the vectors in a frame are in general 
not linearly independent. A frame is called tight if and only if 
the frame bounds A and B are equal. A tight frame is called 
isometric (or uniform) if and only if each vector in F has unit 
norm. 

Theorem 2: Let F be a finite nonempty subset of vectors 
in C d . The following statements about F are equivalent: 

1) F is a uniform tight frame; 

2) F is a WBE-sequence set; 

3) F is a 1-design in C^" 1 . 

Proof: The frame constants of a uniform tight frame F 
in C d are given by A = B = \F\/d, see for example Property 
2.3 in [8]. Therefore, F satisfies equation (0 of Theorem ^ 
for k = 1. The equivalence of the three statements follow now 
from Theorem ^ ■ 
Corollary 1: Any 1-design in CP d_1 is obtained by pro- 
jecting an orthogonal basis from a higher-dimensional space 
(where all basis vectors have the same norm). 



VI. Equivalence of MUBs and 2-Designs 

We need a few more notations before we state our main 
results. If B is a subset of CS d_1 , then the set A = 
{ | (a: | y) | 2 : x,y € B,x ^ y} is called the "angle" set of B. 
For an element x in B and an "angle" a G A, we define the 
subdegree d a (x) as d a (x) = \{y £ B: \(x\y)\ 2 = a}\. If the 
subdegree d a of an a 6 A is independent of x, then B is called 
a regular scheme. Note that the union of mutually orthogonal 
bases of C d is a regular scheme with angle set {0, l/d}. 

Theorem 3: The union X of d + 1 mutually unbiased bases 
in C d forms a 2-design in C5 d_1 with angle set {0, l/d} and 
d{d + 1) elements. 

Proof: We verify that X attains the Welch bound in 
equation Q with equality for < k < 2. The statement then 
follows from Theorem [2 Indeed, this is obvious for k = 0. 
We note that \X\ = d(d + 1). 

If we evaluate the left hand side of the Welch bound for X, 
then we obtain 



x,y£X 



and this coincides with 
Similarly, for k = 2, 

1 



i+i-r 



d(d+l) , s ,1 

- , ' 1+ d-1 0+d 2 - 
d 2 (d+l) 2 V rf 

1 

d 

= l/d; so, X is a 1 -design 



d 2 (d+1) 2 

x,y£X 



d{d+i) 
d 2 (d+i) 2 

2 



l+(d-l)0+cf — 
d- 2 



d(d + 1) 



and this coincides with 



( d+ r) 



= 2/(d(d+l)). ■ 
Theorem 4: A 2-design S in complex projective space 
CS" 4-1 with angle set {0, l/d} and \B\ = d(d + 1) elements 
is the union of d + 1 mutually unbiased bases. 

Proof: A complex projective 2-design with s = 
|{0, l/d} | = 2 satisfies 2 > s — 1, hence is a regular scheme 
[10]. For a = l/d, any x £ B has subdegree d a (x) = d 2 by 
Theorem 2.5 of [10]. It follows that x is orthogonal to d — 1 
elements. 

Let B x = {x} U {z e B: (x\z) = 0}. We claim that B x is 
an orthonormal basis of C d . We may assume that the vectors 
in B are normalized to unit norm. Thus, it suffices to show that 
B x = By for each y G B x . For x = y this is trivial. We know 
that x and y are contained in both B x and B y . Therefore, it 
suffices to show that the intersection set 

I(x, y) = {z€B: (x\z) = 0, (y\z) = 0} = B x n B y - {x, y} 

contains d — 2 elements. 

The number of elements in I(x, y) does not depend on 
x, y for a ^-design with t > 2s — 2, see [11]. Specializing 
Theorem 5.2 in [11] to the case at hand shows that 

l 

\I(x,y)\ =d 2 E <rU°i-M d + VtoM - l - 

ij=0 



We can now evaluate the intersection polynomials gij(0) using 
[11, Theorem 5.3] and obtain that \I(x,y)\ = d — 2. 

Hence we can conclude that each set B x forms an orthonor- 
mal basis of C d . The sets B x partition B. If B x ^ B z , then 
the bases are by construction mutually unbiased. ■ 

Zauner conjectures that if the dimension d is not a prime 
power, then a 2-design with angle set {0, l/d} cannot have 
d(d+l) elements [25]. His conjecture can now be reformulated 
in terms of mutually unbiased bases, which then states that 
N(d) < d + 1 for non-prime power d. If Zauner' s conjecture 
is true, then this would explain the particular role of the finite 
field construction by Wootters and Fields [24]. 

Remark 2: Theorem 3 was obtained earlier by Zauner as 
part of a more general result on combinatorial quantum designs 
using a different terminology, see [25, Theorem 2.19]. The 
converse direction, our Theorem 4, appears to be new. 

VII. SIC-POVMS AND 2-DESIGNS 

Finally, to demonstrate the versatility of TheoremHwe also 
show that another system of vectors used in quantum informa- 
tion theory corresponds to complex projective 2-designs. So- 
called symmetric informationally complete positive operator- 
valued measures (SIC-POVMs) are systems of d 2 vectors 
in C d which have constant inner product, i.e., |(?j,w)| 2 = 
l/(d+ 1) for all v, W in the set. Like in case of MUBs it 
is a challenging task to construct SIC-POVMs — indeed here 
solutions are known only for a finite number of dimensions 
[9], [16]. In [16] it was shown that SIC-POVMs actually form 
complex projective 2-designs. The following theorem gives a 
new and simple proof of this result. 

Theorem 5 (SIC-POVMs are 2-designs [16]): Let X be a 
SIC-POVM X in dimension d. Then X forms a 2-design in 
CS"*" 1 with angle set {l/(d+ 1)} and d 2 elements. 

Proof: Again, we only have to verify that the set X of 
vectors attains the Welch bound with equality for < k < 2. 
The statement then follows from Theorem ^ Indeed, this is 
obvious for k = 0. We note that here \X\ = d 2 . Evaluating 
the left hand side of the Welch bound for X, then we obtain 



x,y£X 



and this coincides with 
Similarly, for k = 2, 



i< 1+ (' i - 1 » = 5 



r;- 1 ) 



l/d; so, X is a 1 -design. 



x,yeX 



1 

d 2 



1 



1 



d(d+ 1) 



and this coincides with ( d+2 L ) = 2/(d(d + 1)). ■ 
Remark 3: Zauner pointed out to us that the previous the- 
orem can also be obtained in the language of combinatorial 
quantum designs by combining Theorems 2.29 and 2.30 in his 
dissertation [25]. 



VIII. Conclusion 

We have shown that the seemingly unrelated concepts of 
MUBs on the one hand and complex projective 2-deigns on 
the other are actually the same objects. This was anticipated 
in a paper by Barnum [2] in which it was shown that the 
union of the (d+1) bases of a particular system of MUBs 
forms a complex projective 2-design. In the present paper we 
have generalized this to arbitrary MUBs and have also shown 
the other direction, i. e., any 2-design in dimension d which 
consists of d 2 + d elements and has angle set {0, 1/d} can be 
partitioned into d+1 sets of MUBs. We have also shown that 
these sets meet the Welsh bounds for k = 0, 1, 2 with equality. 
Hence, the present paper can also be seen as a generalization 
of the results of [21] in which the corresponding statement 
over the real numbers was shown. Finally, we would like to 
mention that Zauner [25] conjectures that affine 2-designs do 
not exist in dimensions d having two distinct prime factors. 
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